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Abstract 

In this paper we investigate the effects that an anomalous acceleration as that experienced 
by the Pioneer spacecraft after they passed the 20 AU threshold would induce on the orbital 
motions of the Solar System planets placed at heliocentric distances of 20 AU or larger as 
Uranus, Neptune and Pluto. It turns out that such an acceleration, with a magnitude of 
8.74 X 10~^° m s~^, would affect their orbits with secular and short-period signals large 
enough to be detected according to the latest published results by E.V. Pitjeva, even by 
considering errors up to 30 times larger than those released. The absence of such anomalous 
signatures in the latest data rules out the possibility that in the region 20-40 AU of the 
Solar System an anomalous force field inducing a constant and radial acceleration with those 
characteristics affects the motion of the major planets. 
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1 Introduction 



The Pioneer 10/11 missions were the first spacecraft to explore the outer Solar System. Pioneer 
10 was launched in March 1972; Pioneer 11 followed it in April 1973. After the encounters with 
Jupiter and Saturn, both Pioneer spacecraft continued to explore the outer regions of the Solar 
System following nearly ecliptical hyperbolic paths in opposite directions. The last successful 
communication with Pioneer 10 occurred in April 2002, at nearly 70 AU from our planet, while 
the last scientific observations were returned by Pioneer 11 in September 1995 when it was at 
approximately 40 AU from the Earth. The last attempt to communicate with Pioneer 10 was 
made on March 2006, but without success ^j. The so-called Pioneer anomaly consists in an 
anomalous constant and uniform acceleration directed towards the Sun of (8.74 it 1.33) x 10"^'^ 
m s~^ found in the data of both spacecraft from the moment when they passed the threshold 
of 20 AU [21 El 15 • Various attempts to explain this feature in terms of known physical and 
engineering effects have been done, but without success; the nature of the Pioneer anomaly is. 
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at present, still unexplained. It is not clear if the anomaly is due to some internal systematic, 
as suggested, e.g., in jSj, or has an external origin. Dedicated space missions have recently been 
proposed [01 El , while the authors of [HI put forth two non-dedicated scenarios involving either 
a planetary exploration mission to the outer Solar System or a piggy-backed micro satellite to 
be launched toward Saturn or Jupiter. 

In this paper we propose to analyze some aspects of the problem of the Pioneer anomaly. In 
particular, we will investigate the possibility that an external, unknown constant and uniform 
force field inducing an acceleration of 8.74 x 10~^° m s~^ on a test particle is present in the outer 
regions of the Solar System within 20-40 AU. The outline of our approach is the following. 

Let us consider a constant and uniform acceleration radially directed towards the Sun and 
of very small magnitude so that it can be treated with the usual perturbative approaches; we 
will deal with it from a purely phenomenological point of view, without making any speculation 
about its origin. For a review of many proposed mechanisms involving possible 'new' physics, 
see |2j and references therein. What is its impact on the orbital motion of the planets of the 
Solar System whose semimajor axis is larger than 20 AU like Uranus, Neptune and Pluto? it 
should be noted that, at present there are no data supporting a more gradual onset of such 
an anomalous effect, and in the framework of this simplified model such planets, due to their 
relatively small eccentricities (euranus = 

0.047, CNeptune — 0.008 and epiuto — 0.248), would be 
acted upon by the anomalous acceleration in all portions of their orbits. 

An analogous approach was recently followed by the authors of JHl- Contrary to us, they 
discard the possibility of using the major planets of the outer Solar System. The focus of ^01 
is on some selected minor bodies characterized by large semimajor axes and eccentricities like 
some unusual minor bodies, Trans-Neptunian Objects (TNOs) and Centaurs. However, such 
objects are not constantly under the action of the hypothetical anomalous acceleration under 
investigation because of their highly eccentric orbits (0.620 < e < 0.947). According to the 
authors of JOl) a relatively low-cost future astrometric campaign should resolve the problem 
within the next twenty years. A different position about the outer planets can be found in [SI. 
Another recent paper in which a comparison between planetary observations and a theoretically 
derived quantity related to the Pioneer anomaly is attempted yielding similar conclusions is . 

In our paper we first work out, both analytically and numerically, the orbital perturbations 
induced by a disturbing acceleration like the Pioneer one on the Keplerian orbital elements 
of a planet. We find that the semimajor axis, the eccentricity, the perihelion and the mean 
anomaly of a test particle are affected by short-periods signatures and the perihelion and the 
mean anomaly undergo also long-period, secular effects. Then, numerical simulations of the 
temporal evolutions of the directly observable right ascensions and declinations of the planets 
induced by a Pioneer-like acceleration are performed. Contrary to the position of the authors 
of ^01 and in according to the conclusions of |S1 , it turns out that the predicted anomalous 
effects of this kind would be sufficiently large to be well detected with the present-day level of 
accuracy in planetary orbit determination jj^, even by considering errors considerably larger 
than those released. Their absence tell us that in the outer region of the Solar System within 
20-40 AU there is no any unknown force field, with the characteristics of that experienced by 
the Pioneer spacecraft, which affects the motion of the major planets. 
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2 The Gauss rate equations 



The Gauss equations for the variations of the semimajor axis a, the eccentricity e, the inchnation 
i, the longitude of the ascending node fi, the argument of pericentre uj and the mean anomaly 
of a test particle of mass m in the gravitational field of a body M can be derived quite 
generally from 

r + ^r = A, (1) 

where fi = G{M + m) and r = — tm- Eq. holds for every perturbing acceleration A, 
whatever its cause or size. The variations of the elements are, for an entirely radial acceleration 



^ = '^'^ Arsinf, (2) 
dc '\J 1 — 6^ 

— = Ar sin /, (3) 

dt na 



duj vT" 



e 



2 



A cos /, (6) 



at na dt 

in which n = 27r/P is the mean motion, P is the test particle's orbital period, / is the true 
anomaly counted from the pericentre, p = a(l — e^) is the semilactus rectum of the Keplerian 
ellipse and A^. is the projection of the perturbing acceleration A on the radial direction of the 
co-moving frame {f,t,n}. As can be noted from eqs. the inclination and the node are 

not perturbed by such a perturbing acceleration. 



3 The orbital effects induced by a constant radial acceleration 



3.1 Analytical calculation 



The ratio of the Pioneer anomalous acceleration to the Newtonian monopole term ranges from 
5 X 10~^ to 2 X 10~^ in the region 20-40 AU, so that eqs. ©-Q can be treated in the perturbative 
way. A treatment of the effects of a constant radial acceleration of arbitrary magnitude in the 
two-body problem can be found in |151 116j . 

The right-hand-sides of eqs. (El)-© have to be evaluated on the unperturbed Keplerian 
ellipse. The secular effects can be obtained by averaging over one orbital period the right-hand- 
sides of eqs. ©-(HI)- 
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The integration of eqs. (dJ-Q is easier if the eccentric anomaly E, defined asAi = E—e sin E 
is used instead of the true anomaly /. The most important relations in terms of E are 

(1 — ecosE') 



dt 
r 

cos/ 
sin/ 



n 



-dE, 



a(l — ecosE), 
(cos E — e) 
1 — e cos E ' 
"v/l — s'mE 
1 — e cos E 



(8) 
(9) 

(10) 
(11) 



By using eqs. (|5))- H11() the shifts of the Keplerian orbital elements over a time span shorter 
than one orbital period, i.e. E — Eq < 2-k, become 



Aa 



Ae 



Aw 



A7W 



2eAr 

— TT- cos^ 
(1 - e^)Ar 



E 
Eo 



COS^ 



Eo 



2 - ^^^0 

-.^ 



n^ae 



Ar 



Eo 
2 sin 2^ 



nAt+ ^-{-SC-e ^ 



+ 



(l-e^)^ 



4e + 



sin^ 



(12) 
(13) 
(14) 
(15) 



Eo 



Prom eqs. ()12() - (|13p it can straightforwardly be noted that the shifts over one full orbital 
revolution, i.e. from Eq to Eq + 27r, vanish for the semimajor axis and the eccentricity, so that 
there are no net secular orbital effects for such Keplerian orbital elements. It is not so for the 
pericentre and the mean anomaly whose secular rates are, from eqs. ()14() - ()15|) 



duj 

m 

dM \ 
dt I 



ArVY 



na 



na 



(16) 
(17) 



The planets of the Solar System show moderate eccentricities (and inclinations). A widely 
used orbital element for such kind of orbits is the mean longitude A = + oj + Ai. Prom eqs. 
(dUl-Jlll) we have, to order 0{e^) 



AA ~ nAt ■ 
The secular rate is, thus 



dX 
H 



1 + — K + -r sin 2^ - -e sin ^ 
4/4 4 



Eo 



~ n 



na 



1 + 



(18) 



(19) 
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3.2 Application to Uranus, Neptune and Pluto 



Since the Pioneer anomaly began to manifest itself at the 20 AU threshold, we will apply such 
analytical results to Uranus, Neptune and Pluto whose semimajor axes are about 19 AU, 30 AU 
and 40 AU, respectively. 

The secular effects which would be induced on the perihelia and the mean longitudes of 
Uranus, Neptune and Pluto by a constant radial acceleration of the same magnitude as that act- 
ing upon Pioneer are listed in Tabled They amount to tens and hundreds of arcseconds/century 
(" cy~^). In regard to their detection, it must be noted that the orbital periods of the investi- 
gated planets, which are 84 years for Uranus, 164 years for Neptune and 248 years for Pluto, are 
comparable or larger than the time span in which accurate observations were collected. Thus, 
it is not yet possible to single out orbital effects averaged over one full orbital revolution for 
Neptune and Pluto. 

In Table [21 which reproduces Table 4 obtained by E.V. Pitjeva in by processing almost 
one century of data of various kinds with the latest EPM2004 ephemerides released by the 
Institute of Applied Astronomy of Russian Academy of Sciences, the formal standard deviations 
of the orbital elements of the nine major planets of the Solar System are reported. The optical 
observations used for Uranus, Neptune and Pluto are 23612 and cover 90 years from 1913 
(1914 for Pluto) to 2003 (see Table 3 of [IHI)- It turns out from Table HI that the accuracy in 
determining the perihelia and mean longitudes of Uranus, Neptune and Pluto would be largely 
adequate to see so huge anomalous signatures, even if we consider that realistic errors could 
be ten-thirty times larger. Indeed, from the formal results of Table |2l we have uncertainties of 
0.1", 2" and 0.3" for the perihelia of Uranus, Neptune and Pluto, respectively; if we multiply 
them by a factor 30 we get 3", 60" and 9", respectively while the anomalous perihelion shifts 
over 90 years would amount to 75.1", 93" and 104.6", respectively. In the case of the mean 
longitudes, the formal uncertainties amount to 0.008", 0.035" and 0.079", respectively which 
become 0.2", 1" and 2" if rescaled by a factor 30; the anomalous shifts on A would be 150.6", 
188.5" and 217.7", respectively. Note also that the systematic errors in the Keplerian mean 
motions 5n = (3/2) y GM/a^Sa due to the (formal) uncertainties in the semimajor axis a amount 
to 0.03 " cy"\ 0.1" cy-i and 0.4" cy"^ only: even re-scaled by a factor 30, they remain well 
smaller than the anomalous shifts. The impact of the mismodelling in solar GM, assumed 
6{GM) = 8 X 10^ m^ s~^ ( http:/ / ssd.jpl.nasa.gov/astro_constants.htm l ) , on the mean longitudes 
is even smaller. 

The very long orbital periods of Uranus, Neptune and Pluto make also the short-period sig- 
nals of interest, in principle, because they would resemble polynomial signals over observational 
time spans like those currently available. This could be useful especially for the semimajor axis 
and the eccentricity which do not exhibit secular rates. In Figure ^ Figure |21 Figure |31 and |1] 
we plot the total, i.e. linear (when present) and sinusoidal, integrated shifts of eqs. (|12|) - (|14|1 
and eq. H18|) for the semimajor axis, the eccentricity, the perihelion and the mean longitude, 
respectively over 90 years. In Table |21 we quote the amplitudes A of the short-period signals for 
Uranus, Neptune and Pluto. Anomalous signatures of such size would be detectable as well, but 
nothing like that is mentioned in Pitjeva's work. 
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4 The orbital effects induced by a constant radial acceleration: 
numerical simulations and comparison with data 



The Keplerian orbital elements are not directly measured quantities: they are related to data 
in an indirect way. For the outer planets the true observables are the right ascension a and the 
declination 5. In Figure 2 of jl3j there are the observational residuals of a cos b and b for Uranus, 
Neptune and Pluto from 1913 to 2004. In order to make a direct comparison with them, in this 
Section we numerically compute theoretical Pioneer-No Pioneer (P-NP) residuals for the same 
quantities over the same time span. 

The first step of this process consists in a direct integration of the gravitational equations in 
Cartesian rectangular coordinates. A number of software packages to perform this goal are freely 
available on the Internet; between them there are OrbFit (http: //newton. dm.unipi . it/orbf it/|). 



ORSA (Orbit Reconstruction, Simulation and Analysis, [http : //orsa ."so urcef orge . net/ by 
P. Tricarico), and Mercury IT7^. The last one has been chosen because it has the simplest use 
and allows to introduce without difficulty user-defined perturbing forces, such the Pioneer one. 
It is introduced through the definition of the corresponding acceleration, whose components in 
X, y and z in an heliocentric frame are given. In this implementation this acceleration has been 
considered only for heliocentric distances greater then 15.0 AU, so to include the whole orbit of 
Uranus, totally excluding Saturn. This package uses the following integration methods 

• second order mixed-variable symplectic 

• Bulirsh-Stoer (general) 

• Bulirsh-Stoer (conservative systems) 

• Radau 15th order 

• hybrid (symplectic / Bulirsh-Stoer) 



Since the pros and cons of a method depend to a certain extent on the specific problem to be 
solved, a number of tests have been performed, using most of the proposed methods and various 
initial conditions in order to choose the most suitable approach to our problem, performed form 
JD 2410000.5 to JD 2460000.5, for the nine planets As a result, the 15th order Everhart-Radau 
algorithm has been chosen for the Pioneer problem and initial conditions at JD 2410000.5. The 
output is in XY Z barycentric, with both integration and tabulation steps of 10.0 days. The 
XY Z coordinates of the Earth's center are obtained from ORSA. The right ascensions and 
declinations of the three planets have been obtained in the standard way (see e.g. [IHj)) and the 
residuals in a cos b and b have been plotted in Figures 151101 As already noted. Figure 2 of [T^ 
shows the observational residuals, in of a cos 5 and b for Jupiter, Saturn, Uranus, Neptune 
and Pluto; the scale is ±5 " and no secular or semi-secular trends can be recognized by visual 
inspection. Our numerically produced P-NP residuals of Figures I5I10I can be straightforwardly 
compared to the last three panels (from the bottom) of Figure 2 of [T^ yielding a direct and 
unambiguous confrontation with the observations. As can be noted, the pattern which would 
be induced by a Pioneer-like acceleration on the planetary motions is absent in the data. 
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5 Summary and conclusions 



In this paper we have considered the orbital motions of Uranus, Neptune and Pluto to inves- 
tigate, from a purely phenomenological point of view, if an anomalous force field inducing a 
constant, uniform Sunward acceleration of 8.74 x 10~^'^ m s~^, as that experienced by the Pio- 
neer spacecraft, affects their motions which occur in the 20-40 AU region of the Solar System, 
well within the range in which the Pioneer data show the anomalous signature. We have worked 
out both analytically and numerically the perturbations induced by such a disturbing accel- 
eration on the Keplerian orbital elements and on the directly observable right ascensions and 
declinations of Uranus, Neptune and Pluto. It turns out that there are long-period, secular rates 
on the perihelion and the mean anomaly and short-period effects on the semimajor axis, the 
eccentricity, the perihelion and the mean anomaly which map onto certain particular patterns 
of the right ascensions and declinations. Such anomalous signatures would be large enough to 
be detected according to the present-day accuracy in orbit determination, but there is no trace 
of them in the currently available data. This result restricts the possible causes of the Pioneer 
anomaly to some unknown forces which violate the equivalence principle in a very strange way 
or to some non-gravitational mechanisms peculiar to the spacecraft. 



Acknowledgements 



L.I. gratefully thanks G. L. Page, D. Izzo and J. Katz for having kindly submitted to his attention 
their useful papers. Thanks also to S. Turyshev for interesting discussions. 



References 

[1] Toth, v., and Turyshev, S.G., The Pioneer Anomaly: Seeking an explanation in newly 
recovered data, gr-qc/0603016 , 2006. 

[2] Anderson, J.D.,Laing, P.A., Lau, E.L., Liu, A.S., Nieto, M.M., and Turyshev, S.G., 
Indication, from Pioneer 10/11, Galileo, and Ulysses Data, of an Apparent Anomalous, 
Weak, Long-Range Acceleration, Phys. Rev. Lett, 81, 2858-2861, 1998. 

[3] Anderson, J.D.,Laing, P.A., Lau, E.L., Liu, A.S., Nieto, M.M., and Turyshev, S.G., 
Study of the anomalous acceleration of Pioneer 10 and 11, Phys. Rev. D, 65, 082004, 
2002. 

[4] Turyshev, S.G., Toth, V.T., Kellogg, L.R., Lau, E.L., and Lee, K.J., The Study of 
the Pioneer Anomaly: New Data and Objectives for New Investigation, jgr-qc/051212l[ 
2005. 

[5] Katz, J.I., Comment on Indication, from Pioneer 10/11, Galileo, and Ulysses Data, of 
an Apparent Anomalous, Weak, Long-Range Acceleration, Phys. Rev. Lett., 83, 1892, 
1999. 

[6] Nieto, M.M., and Turyshev, S.G., Finding the Origin of the Pioneer Anomaly, Class. 
Quantum Grav., 21, 4005-4024, 2004. 



7 



[7] Penanen, K., and Chui, T., A Novel Two-Step Laser Ranging Technique for a Precision 
Test of the Theory of Gravity, Nud. Phys. Suppl. Proc, 134, 211-213, 2004. 

[8] Izzo, D., and Rathke, A., Options for a non-dedicated test of the Pioneer anomaly, 
|astro-ph/05Q4634| 2005. 

[9] Dittus, H., Turyshev, S.G., Lammerzahl, C, Theil, S., Forstner, R., Johann, U., Ertmer, 
W., Rasel, E., Dachwald, B. Seboldt, W., Hehl, F.W., Kiefer, C, Blome, H.-J., Kunz, J., 
Giulini, D., Bingham, R., Kent, B., Summer, T.J., Bertolami, O., Paramos, J., Rosales, 
J.L., Cristophe, B., Foulon, B., Touboul, P., Bouyer, P., Reynaud, S., Brillet, A., Bondu, 
F., Samain, E., de Matos, C.J., Erd, C., Grenouilleau, J.C., Izzo, D., Rathke, A., 
Anderson, J.D., Asmar, S.W., Lau, E.E., Nieto, M.M., and Mashhoon, B., A mission to 
explore the Pioneer anomaly, invited talk given at the 2005 ESLAB Symposium Trends 
in Space Science and Cosmic Vision 2020, 19-21 April 2005, ESTEC, Noordwijk, The 
Netherlands, gr-qc/0506139 , 2005. 

[10] Page, G., Dixon, D.S., and Wallin, J.F., Can Minor Planets be Used to Assess Gravity 
in the Outer Solar System?, astro-ph/0504367j 2005; American Astronomical Society 
Meeting 207, # 154.02, 2005. 

[11] Tangen, K., Could the Pioneer anomaly have a gravitational origin? gr-qc/0602089j 
2006. 

[12] Standish, E.M., An approximation to the errors in the planetary ephemerides of the 
Astronomical Almanac, Astron. Astrophys., 417, 1165-1171, 2004. 

[13] Pitjeva, E.V., High-Precision Ephemerides of Planets-EPM and Determination of Some 
Astronomical Constants, Sol. Sys. Res., 39, 176-186, 2005. 

[14] Milani, A., Nobili, A.M, and Farinella, P., Non-gravitational perturbations and satellite 
geodesy, (Adam Hilger, Bristol, 1987). 

[15] Boltz, F.W., Orbital Motion Under Continuous Radial Thrust, J. of Guidance, Control, 
and Dynamics, 14, 667-670, 1991. 

[16] Prussing, J.E., and Coverstone-Carroll, V., Constant Radial Thrust Acceleration Re- 
dux, J. of Guidance, Control, and Dynamics, 21, 516-518, 1998. 

[17] Chambers, J.E., A hybrid symplectic integrator that permits close encounters between 
massive bodies, Mon. Not. Roy. Astron. Soc, 304, 793-799, 1999. 

[18] Meeus, J., Astronomical Algorithms, (Willman-Bell, 1999). 



8 



Table 1: Secular effects, in arcseconds/century (" cy induced by a constant radial acceleration 
of magnitude 8.74 x 10^^'^ m s^^ on the perihelion and mean longitude of Uranus, Neptune and 



Pluto. 



Planet 


VJ 


A 


Uranus 


-83.5 


167.401 


Neptune 


-104 


209.48 


Pluto 


-116.2 


241.89 



Table 2: Formal standard statistical errors in the non-singular planetary orbital elements, from 
Table 4 of ^Sl- The latest EPM2004 ephemerides have been used. The realistic errors may be 
up to one order of magnitude larger. The units for the angular parameters are milliarcseconds 



(mas). 


Planet 


a (m) 


sinicosri (mas) 


sin i sin (mas) 


e cos w (mas) 


esinru (mas) 


A (mas) 


Mercury 


0.105 


1.654 


1.525 


0.123 


0.099 


0.375 


Venus 


0.329 


0.567 


0.567 


0.041 


0.043 


0.187 


Earth 


0.146 






0.001 


0.001 




Mars 


0.657 


0.003 


0.004 


0.001 


0.001 


0.003 


Jupiter 


639 


2.410 


2.207 


1.280 


1.170 


1.109 


Saturn 


4222 


3.237 


4.085 


3.858 


2.975 


3.474 


Uranus 


38484 


4.072 


6.143 


4.896 


3.361 


8.818 


Neptune 


478532 


4.214 


8.600 


14.066 


18.687 


35.163 


Pluto 


3463309 


6.899 


14.940 


82.888 


36.700 


79.089 



Table 3: Amplitudes A of the integrated short-period effects induced by a constant radial 
acceleration with magnitude 8.74 x 10"^'^ m s~^ on the orbital elements of Uranus, Neptune and 
Pluto. Note that the values of ^aa and the present-day errors in A fully justify the approximation 
used to work out eq. (fTH)) . 



Planet 


^Aa (m) 


^Ae 


^A^ (") 


^AA (") 


Uranus 


1.470 X 10^ 


1 X 10~^ 


237.1 


1.848 


Neptune 


1.03 X 10^ 


2.6 X 10"^ 


3201 


0.82 


Pluto 


6.7 X 10* 


4 X 10"^ 


184.5 


41.64 
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Figure 1: Semimajor axis shifts Aa, in m, for Uranus (a = 19.19 AU), Neptune (a = 30.06 
and Pluto (a = 39.48 AU) over T = 90 years. 
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Figure 2: Eccentricity shifts Ae for Uranus (a = 19.19 AU), Neptune (a = 30.06 AU) and Pluto 
(a = 39.48 AU) over T = 90 years. 
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Figure 3: Longitude of perihelion shifts Aro, in for Uranus (a = 19.19 AU), Neptune 
(a = 30.06 AU) and Pluto (a = 39.48 AU) over T = 90 years. 
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Figure 4: Mean longitude shifts AA, in for Uranus (a = 19.19 AU), Neptune (a = 30.06 AU) 
and Pluto (a = 39.48 AU) over T = 90 years. 
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Figure 5: Numerically propagated P-NP residuals for a cos (5 of Uranus from 1913 to 2004. 
Compare it with the third panel from the bottom of Figure 2 of |13j . 
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Figure 6: Numerically propagated P-NP residuals for 6 of Uranus from 1913 to 2004. Compare 
it with the third panel from the bottom of Figure 2 of |13j . 
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Figure 7: Numerically propagated P-NP residuals for a cos 5 of Neptune from 1913 to 2004. 
Compare it with the second panel from the bottom of Figure 2 of |13j . 
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Figure 8: Numerically propagated P-NP residuals for 5 of Neptune from 1913 to 2004. Compare 
it with the second panel from the bottom of Figure 2 of ^Hj • 
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Figure 9: Numerically propagated P-NP residuals for a cos 5 of Pluto from 1913 to 2004. 
Compare it with the first panel from the bottom of Figure 2 of |13j . 
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Figure 10: Numerically propagated P-NP residuals for 5 of Pluto from 1913 to 2004. Compare 
it with the first panel from the bottom of Figure 2 of jl,3) . 
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